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We propose a scheme that can realize a class of positive-operator-valued measures (POVMs) 
by performing a sequence of projective measurements on the original system, in the sense that 
for an arbitrary input state the probability distribution of the measurement outcomes is faithfully 
reproduced. A necessary and sufficient condition for a POVM to be realizable in this way is also 
derived. In contrast to the canonical approach provided by Neumark's theorem, our method has 
the advantage of requiring no auxiliary system. Moreover, an arbitrary POVM can be realized by 
utilizing our protocol on an extended space which is formed by adding only a single extra dimension. 
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Realization of generalized quantum measurements, or 
positive-operator- valued measures (POVMs), is a funda- 
mental problem in quantum mechanics and quantum in- 
formation science. Many tasks, such as quantum state 
discrimination and entanglement transformation, require 
nonorthogonal measurements to success or to achieve the 
optimal efficiency. During recent years much effort has 
been devoted to the implementation of POVMs on vari- 
ous kinds of physical systems 0]. Many of the proposed 
schemes are derived from Neumark's theorem 2], which 
asserts that any POVM can be realized by extending the 
original Hilbert space to a larger space and performing a 
projective measurement on the extended space. In spite 
of its universality, this method has the drawback of need- 
ing a collective operation on the original system and an 
ancillary system, which may be difficult to implement in 
practice. 

In this letter, we provide a scheme that can realize a 
class of POVMs by performing a series of projective mea- 
surements on the original system alone, in the sense that 
for an arbitrary input state the probability distribution 
of the measurement outcomes is faithfully reproduced. A 
necessary and sufficient condition for a POVM to be re- 
alizable in this way is also derived. Moreover, if only a 
single ancillary dimension is introduced, then arbitrary 
POVMs can be realized by applying our protocol on the 
enlarged space. Nevertheless, our method is limited to 
the physical systems which can be repeatedly measured, 
i.e. we have access to the post-measurement states and 
can perform operations on them again. 

Let us begin with the following example. Suppose we 
are given an arbitrary three-dimensional state, i.e. a 
qutrit. Consider the following protocol: 

(1) Pcrform the projective measurement {|0)(0| + 
|1)(1|, |2)(2|} on the initial state. Suppose the outcome 
r = 0/1 or 2 is obtained. 

(2) If r = 2, do nothing. If r = 0/1, perform the pro- 
jective measurement {|</>o)(</ ) oli |0i)(^>i|> I^X'fel} on the 
state after stage (1), where \4> Q ) = 4?(|0) + |1) + |2)), 



M = ^j(|0)+2|l>-3|2», = _^( 5 |0)-4|1)-|2)). 
Suppose the outcome s — 4>q, (f>\ or 4>2 is obtained. 

(3) If r — 0/1 and s = cf>o, then output the final out- 
come 0; if r = 0/1 and s = (f>i, then output the final 
outcome 1; otherwise, output the final outcome 2. 

Denoting the initial state by p, one can verify that 
the probabilities of obtaining the final outcomes 0,1 
and 2 are tr(pEo),tr(pEi),tr(pE2) respectively, where 
E a = §M))(Vo|, E 1 = &|^i)<M E 2 = I- E - E x 
with |^o) = ^(|0) + II)), = ^(|0) + 2|1)). Thus, 
the above protocol actually can be viewed as a 'virtual' 
POVM {Eo, £1, -E2} if only the measurement outcome is 
concerned. 

The above example demonstrated our basic ideas. In- 
stead of introducing an auxiliary system and perform- 
ing a collective projective measurement on the extended 
system, we here repeat projective measurements on the 
original system over and over again. Suppose we first 
perform a projective measurement {Pi, . . . , P m } on the 
initial state and the outcome i is obtained. Then de- 
pending on i we choose another projective measurement 
{Pi ,...,PmL} and perform it on the state after the 
first measurement. Suppose its outcome is i\. Then bas- 
ing on i and ii we construct another projective measure- 
ment {p[ l,ll \ . . . , Pm'^i} ) } and perform it on the state af- 
ter the second measurement. Similarly the protocol goes 
on. When getting every possible outcome of a performed 
measurement, we may either to output a final outcome 
and finish, or to proceed with a new measurement. All 
these should be specified in advance by the protocol and 
should not depend on the input state. If such a protocol 
outputs the final outcomes with the same probability dis- 
tribution as that of a POVM for an arbitrary input state, 
we say that this POVM is realized by this protocol, al- 
though in our method the final outcome is created more 
artificially. It is worth noting that every protocol has a 
concise graphical depiction in the form of a tree structure 
(For detailed definitions of graphs and trees, see Ref.0]). 
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Each leaf node represents a possible outlet of the proto- 
col and has a corresponding final outcome. For the above 
example, its protocol tree is shown in Fig^ 
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FIG. 1: A graphical depiction of the example protocol. The 
numbers 0/1,2 and symbols 0o> 4>i, 02 along the branches rep- 
resent the possible outcomes of the measurements. The num- 
bers 0,1 and 2 beside the leaf nodes indicate the final outcomes 
inferred from the corresponding chain of measurements. 

In general, we can come up with much more complex 
protocols than the one given above. It is easy to ob- 
serve that any such protocol will generate a POVM. Con- 
versely, we may ask whether an arbitrary POVM can 
be realized in this way. To understand the limit of our 
method, suppose that a POVM {Pi, . . . , E m } is realized 
by a protocol which begins with a projective measure- 
ment {Pi, . . . ,P n }. No matter how the protocol works 
specifically in the subsequent steps, it is always true that 
every POVM element E k should be written as the sum 
of some items M'M where 

y[ _ p{i,ii,—M-i) p(i,h) p(i) p. ^ 

is the product of the projection operators in a chain of 
measurements. Then we have M'M = PJlPi for some 
positive operator II. Thus, Ek could be written in the 
form Ek = y)j—i Pi^kiPi for some positive operators flki- 
It then follows that 

n 

P fe =^PP fc p. (2) 

i=l 

As a consequence, for each E k and each Pi, we have 
[Ek,Pi] = 0. So a necessary condition for a POVM to 
be realizable by our approach is that there exists at least 
one projection operator P ^ I such that P commutes 
with all the POVM elements. 

Next we will prove that the above condition is also suf- 
ficient. Suppose a POVM {Pi, E m } and a projection 
operator P ^ I satisfy that [Ei, P] = for i = 1, . . . , m. 
Let E l0 = PE. t P, E a = {I - P)E l (I - P), and suppose 
they have spectral decompositions 

Eia=f24 a \4aMal (3) 



where a = 0, 1, i = 1, . . . , m and di a = rank(Ei a ) . One 
can verify that Ei = E^ + En, X)"=i ^io — P an d 
Z? =1 E a =I-P : " ' 

Before presenting our protocol for realizing this 
POVM, a lemma should be stated first: 

Lemma 1 If a linear operator M , a state \4>) and a num- 
ber A > satisfy M*M — \\<j)){(j)\ > 0, then there exist a 
state \9) and a number fi > X such that 

A/t|0)=V7#). (4) 



Proof. The proof is given in the appendix. □ 
The constructive proof of this lemma gives us a basic 
function which takes M, X and \<p) as input and outputs 
fi, \9) in the Eq.QJ. We write it in the form 

f(M,X,m = ( f i,\9)). (5) 

Our protocol is as follows: 

Stage 1: Perform the projective measurement {P, / — 
P} on the initial state. If the outcome corresponding to 
P is obtained, set a = 0; otherwise, set a = 1. 

Stage 2: 

(2.1) Set i = 1, j = 1. If a = 0, set M\ Q = P; other- 
wise, set = I — P. 

(2.2) Compute the function 

/«>ALl4 t )) = (/4>IO' (6) 

(2.3) Choose a state |^ a ) e ker(M^). (This is always 
possible, and we will prove it later.) 

(2.4) Perform the projective measurement 

{|^io)<V4U _ IV4>(V4l} 011 the current state, 
where 

V Via V Ra 

If the outcome corresponding to \ip 3 ia ) is obtained, then 
output the final outcome i and exit; otherwise, goto stage 
(2.5). 

(2.5) If j < di a , then set 

Mi^ = {i-m{4 a \)ML (8) 

and increase j by 1; otherwise, set 

Mt i+1)a = (I-\i>lMa\)Ml (9) 

increase i by 1, and set j = 1. 

(2.6) If i = m, then output the final outcome m and 
exit; otherwise goto stage (2.2). 

The protocol can be depicted by the tree shown in 
FigEl One can see that its structure is really simple. 

To prove the validity of the protocol, it suffices to con- 
sider pure input states since the probability of obtaining 
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FIG. 2: A graphical depiction of the general protocol. The 
dotted lines stand for omitted branches. For the measurement 
outcome corresponding to the state \ifP ia ) at every iteration of 
stage (2.4), the symbol ipf a is put along the branch; and for 
the other outcome, we put the performed linear operation 
I ~ IV'iaKV'jal a ^ on S the other branch, including Q a = I — 
^{m-iy^i^^m-Va I where a = 0,1. The numbers 1,. . . ,m - 
1 and m beside the leaf nodes indicate the final outcomes 
inferred from the corresponding chain of measurements. 



each measurement outcome is linear in the density ma- 
trix of the input state. Suppose the input state is It 
can be observed that the role of Mf a is to record the total 
operation performed on the input state, which means, no 
matter at what stage of the protocol, the current state 
is always \tp'} = M- a \ip) /\\M- a \ip)\\ and the protocol can 
reach this stage with probability ||M/ a |i/>)|| 2 . 

Since every Mf a is the product of a chain of projection 
operators which start with either P or / — P, we have 
rank{M 3 ia ) < rank(P) or rank(M^ a ) < rank(I — P), 
which implies ker(M^) ^ 0. So at stage (2.3) the state 
|£| a ) can always be found. 

Consider an arbitrary iteration of stage 2 with i < m. 
First, it is necessary to prove that Mf a , X 3 ia , and \(j> 3 ia ) 
satisfy the condition 

^i-^iwj^ ( io ) 

so that the function / can be applied to them at stage 
(2.2). This will be proved later. Now we assume that it 
holds and then by Eq.QJ get 

MMa) = (11) 



By Eq.0, Eq.lEJ and M/ a f |^ a ) = we have 

Mil\il)&UMi = 4j4 a )(4 a \. (12) 

Then it follows that the probability of obtaining the mea- 
surement outcome corresponding to \ipf a ) is 

MaW)\ 2 = (mil mm^m/mmr 
=K a {m«m a \i>)i\\Mm\\ 2 - 

(13) 

Taking into account the prior probability of reaching this 
stage 1 1 Mf a || 2 , the probability of the protocol stopping 
at stage (2.4) with the current values of i, j and a is 

Therefore, the total probability of the protocol yielding 
the final outcome i is 

1 dia 1 

EE4*<v>i4><4^) = E = (mw) 

a=Q j = l a=0 

(14) 

for all i < m. And naturally the probability of yield- 
ing the final outcome m will be 1 — 1 (V'l^lV') — 
(ip\E m \ip) . So this protocol realizes the POVM 
{.Ei, . . . , E m }. 

Now we go back to prove that the condition (|1C)|) is 
always fulfilled when i < m. Actually, if 

d ia m dj/ a 

MilMi = ^\i\4 a ){4 a \+ e Y,^j4o)(4j 

i'=i+lj'=l 

(15) 

holds, then the inequality (flCjf) will be true. 

We will prove Ea. H15(l by induction on the indices 
We consider only the case of a = 0, because the case of 
a = 1 can be dealt similarly. At the beginning, — 
(1,1), Mi, = P. It follows from Eq.® and E£Li E ™ = 
P that 

771 dio 

EE x L\4oMo\=p- as) 
i=i j=i 

So Ea. holds. Now suppose that for some indices 
Ea. lfTof is valid. If j < d i0 , then by Eq.JBJ and 
Ea. (fT2"|l we have 

dio ., ., ., m d il0 r r r 

= E %o\ti Mo\ + E E ^ol^oX^ol, 
j'=j+i i'=«+ij'=i 

(17) 

which implies that Ea. l|15(l is also valid for the next in- 
dices + 1). Similarly, if j = d i0 , then by Eq.@ 
and Ea. (|12|l the validity of Ea. H15|) for the next indices 
+ can be proved. 

To analyze the efficiency of our protocol, we should be 
aware that its basic idea is to individually realize each 
item ^i a \4^a) Wal m Eq-® and contribute it to the cor- 
responding POVM element Ei for alH < m, while leaving 
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the residual probability to E m . Since each item needs ex- 
actly a projective measurement, our protocol performs at 

m— 1 

most max{ £ rank(E ia ) + 1} projective measurements 

in total. Actually we can rearrange the POVM elements 
{Ex, . . . , E rn } to minimize this upper bound. 
Summarizing, we get the following theorem: 

Theorem 1 A POVM {Ei, E m } can be realized by a 
sequence of projective measurements on the original space 
if and only if there exists a projection operator P =/= I 
such that [Ei, P] = for i = 1, . . . , m. 

As an application, we consider the problem of un- 
ambiguous discrimination 0] of mixed quantum states. 
Suppose a state is secretly chosen from two quantum 
states pi and p2 whose supports have nonempty inter- 
section, i.e. supp(pi) n supp(p 2 ) 7^ 0. Choose a state 
\f) e supp{ Pl ) n supp(p 2 ). If a POVM {E ,E 1 ,E 2 } 
can be used to unambiguously distinguish the two states 
(where E\, E 2 correspond to pi, p 2 respectively, and Eq 
leads to no conclusion), the condition 

tr{E lP2 ) = tr(E 2Pl ) = (18) 

should be fulfilled. Then we have supp(p 2 ) C ker(Ei) 
and supp(pi) c ker(E 2 ), which implies \ip) € ker(E\) 
and \ip) € ker(E 2 ). Let P = \if)){ip\. Then one can verify 
that [Ei,P] = for % = 0, 1, 2. Hence this POVM can be 
realized by our approach. 

A surprising consequence of theorem ^ is that when 
allowing sequences of projective measurements, an arbi- 
trary POVM can be realized by introducing only a single 
ancillary dimension, as the following corollary states: 

Corollary 1 An arbitrary POVM on a d-dimensional 
space can be realized by a sequence of projective measure- 
ments on an extended (d + I) -dimensional space. 

To prove this, note that a POVM {E\, . . . , E m } on a 
e?-dimensional Hilbert space TL can be mapped to the 
POVM {Ex,..., E m , (-01} on any (d + l)-dimensional 
space TL' formed by adding an extra basis element \ip) to 
TL. Letting P = one can find that our condition 

holds trivially. So we can realize this POVM by utilizing 
our protocol on the extended space TL' . Note that in this 
situation, one never needs to perform the projective mea- 
surement {P, I — P} at stage 1, since only the outcome 
corresponding to / — P can be obtained. We can directly 
set o = l. The other part of the protocol remains the 
same. 

In conclusion, we present a protocol that can realize 
a class of POVMs by performing a series of projective 
measurements on the original system, in the sense that 
it can simulate the probability distribution of the mea- 
surement outcomes for any input state. A necessary and 
sufficient condition for a POVM to be realizable in this 



way is also derived. Our method requires no auxiliary 
system and thus may be easier to implement in practice 
than the one provided by Neumark's theorem. Moreover, 
arbitrary POVMs can be realized by adopting our pro- 
tocol on an extended space which is formed by introduc- 
ing only a single extra dimension. Our work may help 
with the implementation of generalized quantum mea- 
surements in the tasks where only the measurement out- 
come is concerned such as quantum state estimation and 
discrimination. 

We gratefully thank Michael Hall for suggesting the 
elegant expression of our theorem using the Lie bracket 
[,] and also pointing out corollary 1. This work was 
supported by the Natural Science Foundation of China 
(Grants Nos. 60503001, 60321002 and 60305005). 

Appendix: Here wc prove lemma ^ Suppose M 
has singular value decomposition M = ^2^ =1 Pi\ipi) (4>i\ , 
where pi > 0, (ipi\ipj) = S i3 , (4>i\4>j) = and d = 
rank(M). Then we have M^M has spectral decompo- 
sition AftM = SLpJl&X&l- By MtM - \\<j>)(<fi\ > 
0, we suppose M'M can also be written in the form 

MtM = X\M4>\+EZiQi\Ci)(Ci\ f o r some m > d- 1, 
qi > and Then by theorem 2.6 of Ref.|5j, we 

conclude that there exists a (to + 1) x (to + 1) uni- 
tary matrix U — (uij)ij—x,..., m +x such that yX\(f>) — 
J2i=x u uPi\<l>i), where ||w|| 2 = J2i=x\ u u\ 2 < 1 - Let 
|0) = Sf=i u «l^)> A* = VIMI 2 > X. Then we obtain 

M*\6) = n^ir £ «iiP<l&> = V7#>- D 

i— 1 



* Electronic address: wgm00@mails.tsinghua.edu.en 
' Electronic address: yingmsh@tsinghua.edu.cn 

[1] E. Andersson, Phys. Rev. A 64, 032303 (2001); S. Franke- 
Arnold, E. Andersson, S. M. Barnett and S. Sten- 
holm Phys. Rev. A 63, 052301 (2001); J. Calsamiglia, 
Phys. Rev. A 65, 030301(R) (2002); S. Virmani and 
M. B. Plenio, Phys. Rev. A 67, 062308 (2003); M. Rosko, 
V. Buzek, P. R. Chouha, and M. Hillery, Phys. Rev. A 68, 
062302 (2003); T. Decker, D. Janzing and M. Rotteler, 
quant-ph/0407054 S. E. Ahnert and M. C. Payne, Phys. 
Rev. A 71, 012330 (2005); M. Ziman and V. Buzek, Phys. 
Rev. A 72, 022343 (2005); S. E. Ahnert and M. C. Payne, 
Phys. Rev. A 73, 022333 (2006). 

[2] A. Peres, Quantum Theory: Concepts and Methods 
(Kluwer Academic Publishers, Dordrecht, The Nether- 
lands, 1993). 

[3] T. H. Cormen, C. E. Leiserson, R. L. Rivest and C. Stein, 
Introduction to Algorithms (Second Edition) (MIT Press, 
2001). 

[4] I. D. Ivanovic, Phys. Lett. A 123, 257 (1987); D. Dieks, 
Phys. Lett. A 126, 303 (1988); A. Peres, Phys. Lett. A 
128, 19 (1988). 

[5] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
2000), p.103. 



